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ABSTRACT 

We study the relative ahgnment of mass and light in a sample of 16 massive early-type galaxies at 
z — 0.2 — 0.9 that act as strong gravitational lenses. The sample was identified from deep multiband 
images obtained as part of the Canada France Hawaii Legacy Survey (CFHTLS) as part of the 
Strong Lensing Legacy Survey (SL2S). Higher resolution follow-up imaging is available for a subset 
of 10 systems. We construct gravitational lens models and infer total enclosed mass, elongation, 
and position angle of the mass distribution. By comparison with the observed distribution of light 
we infer that there is a substantial amount of external shear (7cxt) ~ 0.12, arising most likely from 
the environment of the SL2S lenses. In a companion paper (Ruff et al. 2011) we combine these 
measurements with follow-up Keck spectroscopy to study the evolution of the stellar and dark matter 
content of early-type galaxies as a function of cosmic time. 
Subject headings: galaxies: fundamental parameters — gravitational lensing — 



1. INTRODUCTION 

The last two decades have seen the emergence of a 
standard model of galaxy formation. The main ingredi- 
ents, dominating the overall dynamics of the universe, are 
non-relativistic particles that do not interact with light 
or baryons except through gravity (i.e . cold dark mat- 
ter), and the mysteriou s dark energy ( Perlmutter et al. 
1999 Riess et aL]|1998 ), i.e. a term in the stress-energy 
tensor of the universe characterized by negative pressure 
(or equivalently a cosmological constant). 

The familiar standard model particles move and in- 
teract with each o ther within this ske leton of unkown 
particles and fields ( White & Rees 1978 ) . Although stan- 
dard model particles represent only a s mall minority of 
the t otal energy budget of the universe ( Komatsu et al. 
2011 1, they play a crucial role in the formation of galax- 
les, and their constituents stars and planets. In those 
overdense regions of the universe, interactions between 
standard model particles alter their dynamics and spa- 
tial distribution. Indirectly, through gravitational inter- 
actions, this so-called "baryonic physics" in turn modifies 
the properties of the underlying dark matter distribution 



(Blumenthal et al. 1986 Gnedin et al. 2004 



201Up . Understand ing the interplay between baryons and 



Duffy et al. 



dark matter at sub-galactic scales is crucial not only for 
any effort to understand how galaxies form and evolve, 
but also may shed light on the properties of t he dark 
matter itself (e.g., self-interaction cross- section Spergel 
it] |2000| |Loeb"fc Weiner||2011| ). 
observational point of view, it is very chal- 
lenging to measure accurately the relative distribution 
of baryonic and dark matter on sub-galactic scales. By- 
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and-large, traditional methods require the presence of a 
luminous tracer (e.g. hot plasma, cold gas, or stars), 
whose kinematics are then interpreted to reconstruct the 
underlying gravitational potential (e.g . Saglia et al. 19921 
Gerhard et al.|[200T] [Humphrey et al.n2006; ,lranx et aXl 



19941. Those methods have been used to demonstrate 
beyond any reasonable doubt that galaxies live in dark 
matter halos - if general relativity is a good description of 
gravity at those scales. Gravitation al lensing ( especially 



2010|) provides 
on dark 



strong lensing at sub-galactic scales Treu 
an additional and powerful tool to shed ligh 
matter. By exploiting the deflection of light rays from 
background sources it need not rely on the presence of 
luminous tracers in the deflector. Furthermore, gravita- 
tional lensing is only sensitive to the total gravitational 
potential and therefore can provide accurate measure- 
ments of mass and mass distribution independent of its 
dynamical state or nature. 

The combination of strong gravitational lensing with 
other diagnostic tools, such as stellar kinematics (e.g., 
Miralda-Escud e||l995l iNatarajan fc _Kneib||1996| jTreu k 



Koopmans'2002^ 'Koopmans fc Treu'2003t |Treu fc Koop- 



mans 2004; Koopmans c t al. 2006, 2009), weak lensing 



lations synthesis methods (e.g.. 


Grillo ct al. 


2009 Auger 


et al. 120091 IThomas et al.|20111 


Ibrtora et a 


l.|i>01(]l|Treu 


et al.|2010 Auger et al.||2010), etc, is particularly effec- 



each method alone one can give precise answers to a num- 
ber of questions. What is the relative abundance of dark 
and luminous matter in the inner parts of galaxies? Are 
dark matter density profiles universal as those predicted 
by simulations? Are halos as triaxal as predicted by sim- 
ulations? Are there dark satellites around galaxies? How 
much of the dark matter observed in galaxies is baryonic 
(i.e. low mass stars and high-mass stars remnants) and 
how much is non-baryonic? 

Until a few years ago, answers to these questions based 
on strong gravitational lensing where mostly limited by 
the small samples of known lenses. In fact, strong lens- 
ing is a relatively rare phenomenon, and in general only 
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~ 1/1000 massive galaxy will have a background galaxy 
sufficiently well aligned along the line of sight to pro- 
duce multiple-images. However, this situation is chang- 
ing rapidly owing to the dedicated efforts of a number of 



surveys to find samples of strong lenses ( 


Browne et al. 


2003 


IBolton et al.|2006||Marshall et al.|2009||Treu et al. 


2011 


Brownstein et al. 


2012). In the past few years, 



lenghts have delivered well over 200 galaxy-scale strong 
lenses. 



Th e Strong Lens Legacy Survey (SL2S Cabanac et al 
20071 is a dedicated effort to find strong lens systems 



in the Canada-France-Hawaii Telescope Legacy Survey 
(CFHTLS)^ with the goal of answering fundamental 
questions about the distribution of mass and light. It 
consists of two main efforts: the group and cluster-scale 
survey dLimousin et al.||2009[ [More et al.||2"oIT| [Verdugo 



et al. 
with t 



201 ip and this present series that is concerned 



le follow-up and analysis of the galaxy-scale sam- 
ple (Gavazzi et al. 2012, in prep.). The candidates 
are identified using the RingFinder algorithm which de- 
tects compact rings around centers of isolated galaxies 
(< IO^^Mq), and works by focusing on chromatic image 
excesses around early-type lens galaxies (ETG) that are 
indicative of the presence of lensed arcs. For each of a 
sample of pre-selected bright (Iab < 22.5) red galaxies, 
a scaled, PSF-matched version of the i-band cutout im- 
age was subtracted from the g-band image of the same 
system. The rescaling in this operation is performed 
such that the ETG light is efficiently removed, leaving 
only objects with a Spectral Energy Distribution dif- 
ferent from that of the target galaxy. These typically 
blue residuals are then characterized with an object de- 
tector, and analyzed for their position, ellipticity, and 
orientation, and those showing characteristic properties 
of lensed arcs are kept as lens candidates. A sample of 
several hundred good candidates were visually inspected 
and ranked for follow-up HST imaging HST and VLT 
or Keck spectroscopy that are required to confirm their 
actual lensing nature and to allow accurate lens mod- 
eling. The SL2S galaxy-scale sample provides an ideal 
higher-z complement to the SLAGS: the defiectors have 
similar distribution in size and velocity dispersion, but 
they have a median redshift of ~ 0.5 (c.f. Zd ~ 0.2 for 
SLAGS), extending the baseline for evolutionary studies 
back to half the current age of the Universe. 

In this series we construct lens models for confirmed 
SL2S galaxy-scale lenses, and interpret them by them- 
selves and in combination with stellar kinematics and 



stella r population synthesis models (paper II, Ruff et al. 



2011 1. This paper presents lens models for the 16 con 
firmed systems as well as 1 candidate which proved im 
possible to model with simple gravitational potentials 
and was therefore rejected as strong lenses. We use the 
results to discuss relative orientation of mass and light 
and their flattening. The companion Paper II discusses 
the relative abundance of mass and light and the evolu- 
tion of the mass density profile by combining the SL2S, 
SLAGS, and LSD samples. 
This paper is organized as follows. In Section [2] we 

^ See |htt p: //www. cf ht .hawa li . edu/Sclence/CFHLS/| and links 
therein for a compreliensive description 



introduce the GFHT and HST data used in the analy- 
sis. In Section [3] we describe our modeling techniques. 
Section |4] describes the lens models and presents our re- 
sults on the relative alignment of mass and light. A brief 
discussion of the properties of the lensed sources is also 
presented. Section [5] concludes with a brief summary. 

Throughout this paper magnitudes are given in the 
AB system. We assume a concordance cosmology with 
matter and dark energy density fl,n = 0.3, Qa = 0.7, 
and Hubble constant Ho=70 kms^^Mpc^^. 

2. OBSERVATIONS 
2.1. CFHT data 

The GFHTLS consists of two main components of suf- 
ficient depth and image quality to be interesting for lens 
search*. Both are imaged in the u*, g, r, i and z bands 
with the 1 deg^ field-of-view Megacam Gamera. The 
multi epoch Deep survey covers 4 pointing of 1 deg^ each. 
Two different image stacks were produced: D-85 con- 
tains the 85% best seeing images whereas the D-25 only 
includes the 25% best seeing images. For finding lenses 
we only considered the better resolution stacks. In the 
T06 data release used for this study they reach a typi- 
cal depth of u* ~ 26.18, g = 25.96 ~ 25.47, r ~ 25.43, 
i ~ 25.08 and z ~ 24.57 (80% completeness for point 
sources) with typical FWHM point spread functions of 
0'.'75, 0'.'69, 0'.'64, 0'.'62 and 0'.'61, respectively. The Wide 
survey is a single epoch imaging survey, covering some 
171 deg^ in 4 patches of the sky. It reaches a typical 
depth of u* ~ 25.35, g ~ 25.47, r ~ 24.83, i ~ 24.48 and 
z ~ 23.60 (AB mag of 80% completeness limit for point 
sources) with typical FWHM point spread functions of 
0'.'85, 0'.'79, 0'.'71, 0'.'64 and 0'.'68, respectively Because of 
the greater area, the Wide component is our main source 
of lens candidates. 

Around each lens and in each GFHTLS band, we pro- 
duce cutout images 101 pixels {i.e. 18'.'8) wide. Nearby 
stars are used to produce a synthetic PSF model. The 
17 systems analyzed in this study are shown in Figure [l] 
and listed in Table O 

2.2. HST follow-up imaging 

In order to confirm the lensing hypothesis and allow 
for detailed lens modelling, 65 galaxy-scale lens candi- 
dates have been observed with the Hubble Space Tele- 
scope (HST) as snapshot programs during cycles 15, 16, 
17 (GOs 10876, 11289 PI J.-P. Kneib; GO 11588 PI 
Gavazzi). The observations started with the Advanced 
Gamera for Surveys (AGS), then switched to the Wide 
Field and Planetary Gamera 2 (WFPG2) after the failure 
of AGS, and finally turned to the Wide Field Gamera 3 
(WFG3) after Servicing Mission 4. Approximately 50% 
of the lens candidates were confirmed as lenses in this 
way. A more comprehensive description of the efficiency 
of the SL2S lens finding strategy is given by Gavazzi et 
al.(2012, in preparation). 

Ten of the galaxy-scale systems were observed with 
HST early enough to be included in our first Keck spec- 
troscopic follow-up campaign and are the subject of pa- 
per I and II. The remaining systems will be presented 
in a forthcoming paper at the end of the Keck and HST 

See http://t erap lx.iap.fr/cplt/T0006/T0006-doc.pdf I 
| |Goranova et al.|2009[ l 
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TABLE 1 

Summary of HST observations. 



Fig. 1. — Gallery of the 17 SL2S systems studied here (composite rgb color obtained with Megacam g, r and i filters). Prom top left to 
bottom right, the lenses are sorted like in table [2] Images are IS'.'S on a size and oriented North up and East left. 

Exposure times and final pixel scales are given in Table[T] 

3. MODELING METHODOLOGY 

In order to model the light distribution of a galaxy- 
scale gravitational lens, one has to disentangle the con- 
tribution of the foreground deflector and that of the back- 
ground lensed arc-like features. The former is generally 
a red ETG and the former a dimmer, blue, and presum- 
ably star-forming, more distant source. The light dis- 
tribution of an ETG generally ha s a r egular shape suffi- 
ciently well described b y a Sersic d Sersic 1968) or even a 
de Vaucouleurs profile (de Vaucouleurs 1948) with very 
small color gradients and blue background sources are 
also well represented by Exponential profiles. The sep- 
aration of these two components is relatively straight- 
forward for ETGs deflectors and they can generally be 
fitted independently for most applications. However, 
very detailed investigations of the source properties of 
high signal-to-noise data mi ght require a simult aneous 



name 




instrument / filter 


exp. time 

(sec) 


pixel scale 

(arcsec) 


SL2SJ021411 


-040502 


ACS/F606W 


400 


0.05 






ACS/F814W 


800 


0.05 


SL2SJ021737 


-051329 


ACS/F606W 


400 


0.05 






ACS/F814W 


800 


0.05 


SL2SJ022511 


-045433 


WFPC2/F606W 


1200 


0.1 


SL2SJ022610 


-042011 


WFPC2/F606W 


1200 


0.1 


SL2SJ022648 


-040610 


WFPC2/F606W 


1200 


0.1 


SL2SJ023251 


-040823 


WFPC2/F606W 


1200 


0.1 


SL2SJ140123- 


f555705 


WFPC2/F606W 


1200 


0.1 


SL2SJ141137+565119 


WFC3/F475X 


720 


0.04 






WFC3/F600LP 


720 


0.04 


SL2SJ221326 


-000946 


WFPC2/F606W 


1200 


0.1 


SL2SJ221407 


-180712 


WFPC2/F606W 


1200 


0.1 



observing campaigns. Of these ten systems, three were 
observed in two bands, with ACS or WFC3 whereas the 
others were observed only in a single WFPC2 band as 
detailed in Table [T] 

All the WFPC2 data were reduced using standard 
MultiDrizzle^ recipes. The cosmic ray removal worked 
well because the 1200 sec exposure time was split into 
3 exposures. The "drizzling" was performed by preserv- 
ing the CCD frame orientation and pixel scale to avoid 
producing correlated noise. The ACS and WFC3 obser- 
vations consisted of single or double exposures onl y and 



fit of both components {e.g. Marshall et al.|[2"007 ), espe 



cially with ground-based data (Brault et al., m prep) and 
bright arcs. We describe our procedure to fit the fore- 
ground defiector in |3.1[ In Section [3?2| we describe how 
we use the residual lensed images to model the gravita- 
tional potential of the deflector and the intrinsic surface 
brightness of the source. 

3.1. Foreground deflector 



therefore we relied on LA-Cosmic (van Dokkum 2001 1 



on individual exposures for c osmic ray removal, before 
combining them with swarp ( Berlin et al. 2002 ) . Like 
for WFPC2 images are kept in the natural CCD frame. 

® MultiDrizzle is a product of the Space Tel escope Science In- 
stitute, which is operated by AURA for NASA | |Koekemoer et aL] 
|2002[ |. 



We used the versatile galfit software (Peng et al 
2002 1 to perform the subtraction of the foreground de- 
flector as it allows to account for boxy-/diskyness of 
isophotes. It generally yields good image subtraction 
with a Sersic profile (see Fig. [4]) but the recovered Sersic 
indices n and effective radii Ttcff are quite degenerated 
and this often leads to misleading values of n and i?cff- 
We thus used a generic Sersic profile to get as good a 
deflector subtraction as possible and we forced n = 4 to 
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get more robust values of n and Rett- In the same vein, 
the fit of the n — 4 foreground lens also yields a robust 
measurement of ellipticity and orientation of the light 
distribution. These values are reported in table [2] 

Fits are performed in all the available bands of a given 
lens using 18'.'8 on a side cutout images and suitable 
PSF eith er inferred from th e neighboring stars or from 



TinyTim (jKrist et al. 



2011 



10 ; 



in the case of HST data. 

The formal errors on each parameter are generally very 
small due because they only account for statistical errors 
and not modelling errors associated with the somewhat 
over-simplistic assumptions about the light distribution. 
A more realistic estimate of the total uncertainties on the 
recovered ellipticity, orientation and effective radius can 
be made by estimating the filter-to-filter dispersion on 
these parameters when HST data is missing. When HST 
imaging is available, shape parameters characterizing the 
deflector are more robustly measured and we adopt the 
formal errors from galf it. 

3.2. Lensed features and mass modeling 

For lens modehng we used a dedicated code sl_f it de- 
veloped for and teste d on galaxy- scale strong lens es (e.q. 



Gavazzi et al.l|2007| |2008, ,2011i |Tu et al. |2009|). The 
code hts model parameters of simple analytic tensing po- 
tentials. It uses the full surface brightness distribution 
observed in the image plane and attempts to explain it 
with one or more simple analytic light components de- 
scribed by an exponential radial profile with elliptical 



shape (see e.g. [Marshah et al.||2007[ [Bolton et aL||2008 
for similar techniques). 

The lensing potential is assumed to be made of a sin- 
gular isothermal ellipsoid (SIE), centered on the main 
defiector. This is the simplest mass profile that has 
been shown to yield a description of the mas distribution 
of massive early-type galaxies sufficient to derive Ein - 
stein Radius, position angle and ellipticity (Treu 2010) 



The convergence profile of the central mass component 
is given by 

i^ix,y)^^, (1) 

where the scaling parameter b is the Einstein radius 
i?Ein and is related to the velocity dispersion of the de- 
flector through 6/1" = (cr„/186.21 kms^^fDis/D^ and 
= qx^ + /q is the radial coordinate that accounts for 
the ellipsoidal symmetry of the iso-density contours and 
q is the minor-to-major axis ratio. The orientation of the 
major axis PAtot is allowed to vary, although this is not 
explicit in the definition of Equ. Q . We do not assume a 
priori that the orientation of the total mass distribution 
is correlated to that of the observed stellar component. 
This way the total potential can capture (and mix be- 
cause of substantial degeneracies) i nternal and externa l 
quadrupolar terms in the potential ( Keeton et al.|[T997 ). 
We also add an external shear parameter when more 
information is available (see the case of SL2SJ021737- 
051329 described below), sufficient to break the degener- 
acy between external shear and orientation of the galaxy 
potential. 

Like in the previous step, when we fitted the fore- 
ground light distribution we weigh pixels with the im- 



age inverse total variance (including sky, foreground and 
lensed features) . The term relating the observed light 
distribution I{xi) at pixel i and the intrinsic source light 
distribution S{y) is thus: 



X 



E 



[lixi) - S{xi- d (xi|pp) \ps)Y 



(2) 



where = [6,9,PAsie] contains the parameters 

that determine the gravitational potential and Ps — 
[xs, ^s, (/s, PAg, TOs, Tg] is the section of the parameter 
space that codes for the source light distribution, namely 
its center Xs, its ellipticity q^^ its orientation PAg, its 
magnitude rus and half-light radius Tg. 

The optimization of these 9 parameters is performed 
using Monte Carlo Markov Chain techniques. Reported 
model parameters and confidence intervals are taken 
from the 0.50, 0.16 and 0.84 quantiles. However, be- 
cause of the simplicity of both the model potential and 
model source light distribution, we generally end up with 
very small formal errors on the recovered model parame- 
ters th at should be substantially increased. [Bolton et al. 
(20081 estimated that relative errors on i?Ein should be 
about 5%. We thus add this dominant contribution in 
quadrature to the statistical errors in table [2| and do the 
same for the axis ratio gsiE- 

4. RESULTS 

In this section we describe the lens model of each sys- 
tem (§ 4.1). In addition, we study the shape and relative 
orientation of the light and total mass distributions as 
they do are independent of the source and deflector red- 
shift (§ 4.2). Results depending on spectroscopic infor- 
mation (source and deflector redshifts and velocity dis- 
persion) quantity along with a novel method proposed to 
mitigate the nuisance due to the ignorance of the source 
redshift are presented in Paper II. 

4.1. Notes on Individual Lens models 

For each lens in the sample we describe the resulting 
best fit model: 

• SL2SJ021411-040502 is well reproduced by the sin- 
gle SIE potential + Exponential source model. We 
note a small axis ratio gsiE 0.33 this seems 
at odds with the more circular light distribution 

~ 0.89. This can be explained by the presence of 
a neighbor galaxy at the same redshift as Zd about 
8" and PA ~ 65° which enhances the elongation of 
the potential. 

• SL 2SJ021737- 51329 has been studied by in detafl 
by Tu et al. (2009). Our results are in agree- 
ment with theirs. A simple SIE -I- Exponential 
source yields a good fit, even though a better fit 
is achieved by introducing both internal elliptic- 
ity through the SIE potential and external shear 
(consistent with the presence of a nearby group of 
gal axies) . Similar concl usions were recently found 
by Cooray et al. (2011) using additional Near IR 
CAN DELS data. TEesource redshift is = 1.847 
(Paper II). This implies that the source half-light 



radius is i?cff,s 



0"08il^:°?« 



0.69 



+0.24 
-0.13 



kpc for 



http: //www. stsci . edu/hst/observatory/f ocus/TinyTim 



an F606W magnitude = 24.45 ± 0.10. The sys- 
tem is lensing a second source at Zs2 ^ 2.3 that 
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Fig. 2. — Results of the lens modeling analysis. For each row, a system is shown with, from left to right, lensed features with the deflector 
subtracted off with galf it (data), the best fit model prediction (model), the residual (data-model) and the associated source plane light 
distribution. The critical lines are overlaid in the first three panels whereas the caustics lines are shown in the last one. 



we do not consider here for consistency with the 
other systems. However, regardless of whether the 
second source is modelled or not, the good signal- 
to-noise ratio and the favorable image configura- 
tion allows us to break the degeneracy between ex- 
ternal shear and internal cUipticity and we there- 
fore get constraints on both. We find an exter- 
nal shear 7ext = 0.064 ± 0.003 with an orientation 



PAext = 89?7 ± 0?4. We sec that the light is well 
aligned with the external shear and about ^ 33 deg 
misaligned with the SIE component. This may be 
due to tidal effects. 



SL2SJ021902-082934 is another "cusp" configura- 
tion with a marginal candidate counter image on 
the opposite side. The lack of HST imaging im- 
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SJ022610-042011 
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Fig. 3. — Lens models (continued). 

plies that we had to consider the CFHT g band 
image. We can see an additional source westward 
of the arc that has a color similar of the deflector 
and should not be considered as a lensed feature. 
Since we can achieve a good fit of the arc without 
introducing this perturbation we neglect it. 

• SL2SJ022056-063934 is a minor axis cusp config- 
uration. The single source component does not 




REboNSTRudTED SOURCE 





REboNSTRudTED SOURCE 




REboNSTRudTED SOURCE 




perfectly capture the faint tail of the inner arc al- 
though most of the fiux is well recovered, even with 
ground-based resolution. 

• SL2SJ0225 11-045433 is a low redshift bright de- 
fiector with another bright minor axis cups con- 
figuration. The source redshift is Zg = 1.199 (pa- 
per II). We report here the result with a single 
source component of magnitude nis = 24.14 ± 0.04 
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Fig. 4. — Lens models (continued). Note the different plotting range for the first row corresponding to SL2SJ140614+520252. 



in the F606W band and half-light radius i?cfi' s — 
0'.'125 ± (y.'OOS ~ 1.037 ± 0.025 kpc. We note we 
also attempted to account for the fainter exten- 
sion of the furthest arc (on the East of the de- 
flector), with a secondary component. This would 
even lower the residuals on the two multiple images 
without changing the results on the recovered po- 
tential parameters. Accounting for this component 
the source would be ~ 0.36 magnitudes brighter. 



• SL2SJ022610-042011 is a typical large impact pa- 
rameter double configuration implying a substan- 
tial differential magnification of the two multiple 
images. In addition, these two are nearly aligned 
with the center of the deflector. This is consistent 
with the potential being close to circularly sym- 
metric (gsiE ^ 0.92). This system has a known 
source redshift = 1.232. We can thus estimate 
the source half-fight radius i?eff,s = 0'.'160±0'.'011 ~ 
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Fig. 5. — Lens models (continued). Note that the system in the first row, SL2SJ221606-175131, is not a satisfying model of a lens and 
we disqualify it as being an actual gravitational lens. 



1.332 ± 0.092 kpc for an F606W magnitude = 
25.10 ±0.08. 

• SL2SJ022648-040610 is another double configura- 
tion with a more balanced magnification ratio be- 
tween the two images involving a slightly more 
elongated potential. We notice that the deflector 
looks like an edge-on SO galaxy with an elongated 
light distribution. However, the total potential is 
considerably more circular. 

• SL2SJ022648-G90421 is also modeled using Mega- 
cam g band data and shows a minor axis cusp con- 
figuration with little deviation between light and 
SIE orientations. This is the faintest source we re- 
construct with TTig = 27.13 ± 0.14. 

• SL2SJ023251-040823 is a double system with the 
source close to the cusp. 

• SL2SJ140123-I-555705 shows a dim cusp-like arc. 
Even though we cannot identify the counterimage, 
the substantial bending of the arc breaks the degen- 
eracy between shear and intrinsic source ellipticity 
and allows us to measure the Einstein radius with 
good accuracy i?Ein = 1''19 ± O'.'OT, comparable to 
the other cases. 

• SL2SJ140614-I-520252 is another cusp configura- 
tion on a larger scale, as i?Ein = 3'.'02±0"15. Some 
blue light excesses are seen in the g band image 
shown the Fig. |4] and only the main "naked" cusp 
arc is captured by our model. However we can see 
that the other Western component is not multiply 
imaged as the source would fall outside the caustic. 
Thus it does not provide more information and we 
verified that its inclusion does not change our con- 
clusions on the lensing potential. Even by ignoring 
this second singly imaged source component, we 
find that this system involves the brightest source 
of the sample with g band nis = 23.50 ± 0.03. We 



could not find a solution with a lensed source re- 
producing the Southern and Northeastern blue ex- 
cesses. The lens potential is highly elongated with 
gsiE = 0.290 ± 0.015. This is potentially due in 
part to shear from a galaxy about 8'.'3 in the South 
East direction. 

• SL2SJ141137-F565119 is the only system having 
WFC3 imaging presented here. The signal-to-noise 
ratio is good and allows us to get tight constraints 
on the potential and on the source whose redshift 
is found to be Zg = 1.420. The source half-light 
radius in the F475X band is ReS,s = 0'.'0.058 ± 
O'.'OOl ~ 0.490 ± 0.008 kpc and its magnitude is 
TOs = 25.86 ± 0.04. We note that the addition of 
a very compact core would improve the fit further 
as our single source component does not reproduce 
the full complexity of the source light distribution. 
Our best fit model nevertheless captures most of 
the extended Einstein ring. 

• SL2SJ220629-h005728 is modeled using the g band 
Megacam image. However, we see in Fig. [TJthat a 
small red satellite lies on top of the lensed blue fea- 
tures. We thus need to disentangle both contribu- 
tions by simultaneously fitting to the g and i bands 
one lensed source and one foreground unlensed 
source centered in the red satellite. This satel- 
lite also contributes to the lens potential through 
a point mass component of unknown mass. We 
find that this perturbing mass has to be Mport ^ 
2 X 10^ Mq (68% CL). Getting HST imaging data 
would significantly improve the accuracy of the de- 
composition and the constraints on Mpert- The 
Einstein radius estimate is robust with respect to 
the inclusion (or not) of the perturbing potential. 

• SL2SJ221326-000946 is an edge-on disk galaxy in 
which both the potential and the light distribu- 
tion are highly elongated: qsiE = 0.191 ±0.014 and 
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= 0.277 ± 0.007 with a small APA ■ 
misalignment. 



10?4±0?6 



• SL2SJ221407-180712 does not exhibit any counter- 
image nor a strong curvature of the main outer im- 
age. This implies that the model cannot be con- 
strained well and we can only place upper limits 
on the Einstein Radius i?Ein = 0'.'41 ± 0'.'23. The 
results on the potential elongation and orientation 
are thus very loose and we do not consider this sys- 
tem in the statistical analysis of { 4.2 Note that the 
small upper limit on i^Ein is consistent with other 
mass constraints in paper II. 

• SL2SJ221606-175131 has a relatively symmetric 
configuration. Its similarity to a classic quad con- 
figuration lead us to include it in the lens candidate 
sample. However, we could not find any sensible 
lens model able to reproduce in detail the image 
configuration. Therefore we conclude that this sys- 
tem is not a lens, or more precisely, that the blue 
features about 1" from the center of the ETG do 
not originate from a unique background source. We 
speculate tha it could be due to low surface bright- 
ness star formation at the ETG's redshift. We thus 
exclude this system from the statistical analysis be- 
low. 

• SL2SJ221929-001743 does not have HST imag- 
ing and the CFHT data do not allow us to con- 
straint tightly the potential parameters as we can- 
not identify unambiguously a counterimage. In 
addition, the arc does not display a strong cur- 
vature at the resolution of ground-based imaging. 
We thus get weaker constraints on the Einstein 
Radius i?Ein = 0'.'74 ± 0'.'14 than for other sys- 
tems. HST imaging would significantly improve 
the mass model. Given the redshift of the source 
Zs = 1.023, we infer a source half-light radius 



cff.s 



— n"'?»/1+0073 



3.03_o;43 kpc and a g band 
intrinsic magnitude TOs — 24.08 ± 0.21. 

4.2. Alignment of mass and light 

The top panel of Figure [6] compares the ellipticity of 
the stellar component to that of the mass model. For the 
vast majority of the objects the cUipticitics arc tightly 
correlated. This is expected since stellar mass makes up 
a significant fraction of the total mass within the Ein- 
stein Radius. However, there are few interesting out- 
liers. The flattened light distribution of SL2SJ022648- 
040610 is consistent with a disky galaxy living in sig- 
nificantly rounder halo. Perhaps more suprising arc the 
cases of SL2SJ021411-040502 and SL2SJ220629-h005728, 
where the stars are almost round, while the potential is 
significantly flattened. In both cases the source of the 
ellipticity of the mass distribution seems to be external 
shear, associated with a satellite. In general, the relation 
between light and mass ellipticity seems to have signifi- 
cantly more scatter than that found for the SLAGS sam- 
ple, < qsis/q* >— 0.9 5 with scatter 0.48, as opposed to 
0.99 with scatter 0.11 ([Koopma ns et al.|2006 1. As we dis- 
cuss in the next paragraph this is probably dire to a much 
more significant role played by external shear. The bot- 
tom panel of Figure [6] shows the angular offset between 



I 
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Fig. 6. — Top panel: correlation between axis ratio of the 
light and of the mass distribution. Bottom panel: angular off- 
set between the major axis of the mass distribution and that 
of the light as a function of axis ratio of the mass distribution. 
SL2SJ021737— 051329 is not included because this system was 
modeled including external shear. 

light and mass PA as a function of mass ellipticity. Even 
neglecting the points with gsiE > 0.85, where clearly the 
position angle is not very well measured since the po- 
tential is so circular, there is considerable scatter around 
zero. For the entire sample the rms scatter is 25 degrees, 
while if we limit ourselves to qsie < 0.85 the scatter is 
still 18 degrees. Again this is considerable higher than 



the 10 degr ees found for the SLAGS sample (Koopmans 
et al. 20061 in general, and closer to the values found 



for the subset of SLAGS lenses that reside in overdense 
environments ( |Treu et al. 2009 1). A simple calcula tion, 
based on Eq 22 in the paper by |Keeton et al. ( 1997[ ) and 
assuming that the direction of external shear is randomly 
distributed with respect to that of the mass of the galax- 
ies shows that this is consistent with a relatively large 
average external shear. This is illustrated in Figure [7| 
where we show the expected r.m.s. fluctuation of the 
position angles as a function of average external shear. 
Regardless of whether the rounder objects are included 
or not, it seems that external shear of order « 0.12 is 
required on average. This level of external shear is fairly 
common amongst galaxy-scale gravitational l enses (e.g. 



Keeton et al. 1997 



Holder fc Schechter||200 3 [), and it is 



ikely due to the environment of the lenses (fXuger et aL 



2007 Treu et al. 



since massive ear 



2009} |Auger|[2008l ^Wmg et al.||2011^ 



y-type galaxies typically reside at the 
centers of groups. The lower level of external shear in the 
SLAGS sample is most likely due to the smaller size of 
their Einstein radius relative to the characteristic scale of 
the galaxy (half effective radius typically) , and therefore 
the more relevant role played by stellar mass in defin- 
ing the potential within the critical curve. We note also 
that, despite a greater influence of the environment rel- 
ative to SLAGS, the SL2S lenses do not require as large 
an quadrupolar term as JVAS radio source or quasars 



(Keeton et al. 1997). This might be due to the fact 



that extended sources are not as sensitive to the mag- 
nification bias that would boost the fraction of highly 
magnified quads. In turn, the only unambiguous quad 
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Einstein radii is gre ater that for SLAG S: i?Ein/-RGff' — 1-1 

In addition, we 



2010) 



0.06 0.08 

Fig. 7. — The curves represent the expected rms fluctuations 
of the ofl'set in PA as a function of average external shear. The 
horizontal lines represent the values measured for the SL2S sample. 
The solid lines refer to the whole sample, while the dashed lines 
refer to the subsample with mass axis ratio gsiE < 0.85. 

systems here are J141137+565119 and J021737-051329. 
Two out of 16 lenses (excluding J221606— 175131) are 
quads. 

5. SUMMARY 

In this paper we have presented gravitational lens mod- 
els of a pilot sample of 16 galaxy-scale lenses identified as 
part of the SL2S Survey. Hubble Space Telescope imag- 
ing is available for 10/16 systems, while for the others the 
modelling is based on ground based GFHTLS imaging. 
After removing the light from the foreground deflector we 
use the surface brightness distribution of the background 
source to constrain mass models of the deflector galaxy. 
For each system we derive Einstein Radius, position an- 
gle and ellipticity of the mass distribution described as 
a singular isothermal ellipsoid. These parameters are 
used in paper II in combination with spectroscopic infor- 
mation to study the relative distribution of stellar and 
dark matter in early-type galaxies and its evolution with 
cosmic time. In this paper we focus on the relative ori- 
entation and ellipticity of the luminous and total mass 
distribution, which does not require spectroscopic infor- 
mation. We find that ellipticity of mass and light are 
tightly correlated except for a handful of outliers. In 
one case the presence of a disk makes the light distri- 
bution significantly flatter than the overall mass distri- 
bution. In two cases, the presence of a nearby galaxy 
introduces significant external shear in the overall mass 
distribution. In addition, we find that the position an- 
gle of mass and light is on average aligned, albeit with 
r.m.s. scatter of 18-25 degrees, significantly larger than 
what is found for the lower redshift SLAG sample. We 
interpret this scatter as evidence for substantial external 
shear, on average (7oxt) ~ 0.12, resulting from the envi- 
ronment (physically related to the main defiector galaxy 
or along the line of sight) The physical scale of the SL2S 



vs 0.5 for SLAGS ( Auger et al.] 

find only two quad lenses in this pilot sample, suggesting 
that magnification bias is less effective in boosting the 
statistics of extended lensed sources as compared with 
samples comprising point sources. Another difference be- 
tween SL2S and other samples is that SL2S is a deflector 
selected sample, as opposed to source selected samples 
like GLASS. 

In paper II we combine our lensing information with 
stellar kinematics to infer the cosmic evolution of the 
mass density profile of massive galaxies. Given the sam- 
ple size in paper II, the role of environment and exter- 
nal shear cannot be explored further, neither in terms 
of astrophysical signal, nor in terms of systematic uncer- 
tainty on the mass density slope. However, the follow-up 
imaging and spectroscopy of SL2S candidates is still on- 
going and we plan to investigate these issues with a much 
larger sampe in future papers of this series. The paral- 
lel effort of measuring the redshift of sources with Keck 
and XShooter on VLT will also allow us to explore the 
population of lensed sources with greater detail. 
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